Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this work, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.     

Imposing accurate wall boundary conditions in corrective‐matrix‐based moving particle semi‐implicit method for free surface flow

Corrective matrix that is derived to restore consistency of discretization schemes can significantly enhance accuracy for the inside particles in the Moving Particle Semi‐implicit method. In this situation, the error due to free surface and wall boundaries becomes dominant. Based on the recent study on Neumann boundary condition (Matsunaga et al, CMAME, 2020), the corrective matrix schemes in MPS are generalized to straightforwardly and accurately impose Neumann boundary condition. However, the new schemes can still easily trigger instability at free surface because of the biased error caused by the incomplete/biased neighbor support. Therefore, the existing stable schemes based on virtual particles and conservative gradient models are applied to free surface and nearby particles to produce a stable transitional layer at free surface. The new corrective matrix schemes are only applied to the particles under the stable transitional layer for improving the wall boundary conditions. Three numerical examples of free surface flows demonstrate that the proposed method can help to reduce the pressure/velocity fluctuations and hence enhance accuracy further.     

石墨基液流电池复合双极板

液流电池是一种安全性高、使用寿命长、可扩展的大规模储能系统,可以协助电网调峰储能,提高能源利用率,发展前景广阔。双极板是液流电池的重要组成部分。功能上起到了分隔、串联电池、传导电流、为电堆提供结构支撑等作用。从成本构成角度看,双极板的价格占电堆成本的比重也较大。开发高性能、低成本的双极板对加快液流电池的商业化应用具有重要意义,也是目前业界的迫切需求。虽然文献上报道了许多针对液流电池双极板开发的工作,但是目前高性能、低成本的液流电池双极板产品仍无法充分满足市场需求。本文着重介绍了石墨基复合双极板的研究现状,介绍了材料选择、工艺流程对关键性能的影响,对相关工作进行了评述,并为液流电池双极板的开发提出了建议。     

Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

Cytotoxic dimeric half-sandwich Ru(II), Os(II) and Ir(III) complexes containing the 4,4′-biphenyl-based bridging ligands

A series of dinuclear half-sandwich Ru(II), Os(II) and Ir(III) complexes [Ru₂(μ-Lⁿ)(η⁶-pcym)₂Cl₂](PF₆)₂ ( 1 , 4 ), [Os₂(μ-Lⁿ)(η⁶-pcym)₂Cl₂](PF₆)₂ ( 2 , 5 ) and [Ir₂(μ-Lⁿ)(η⁵-Cp*)₂Cl₂](PF₆)₂ ( 3 , 6 ), based on 4,4′-biphenyl-based bridging Schiff base ligands N,N′-(biphenyl-4,4′-diyldimethylidyne)bis-2-(pyridin-2-yl)methanamine (L¹; for 1 – 3 ) and N,N′-(biphenyl-4,4′-diyldimethylidyne)bis-2-(pyridin-2-yl)ethanamine (L²; for 4 – 6 ) is reported; pcym = 1-methyl-4-(propan-2-yl)benzene, Cp* = pentamethylcyclopentadienyl. The complexes were characterized by relevant analytical techniques (i.e. elemental analysis, FT-IR, NMR, ESI-MS), and their in vitro cytotoxicity was assessed at six cancerous and two non-cancerous (healthy) human cell lines. Overall, complexes 4 – 6 , containing the L² bridging ligand, revealed higher cytotoxicity as compared with 1 – 3 and, thus, they were studied in greater detail. The best-performing complex 6 exceeded at least twice the in vitro cytotoxicity of cisplatin and showed high selectivity towards the cancer cells over the normal ones, including the primary culture of human hepatocytes. In contrast to cisplatin, complexes 4 – 6 did not induce the cell cycle modification of the treated A2780 human ovarian carcinoma cells (studied by flow cytometry and Western blot analysis). High levels of superoxide anion were induced by complexes 4 – 6 at the A2780 cells. The levels of activated forms of Caspase-3 and Caspase-8 at the A2780 cells treated by Ru(II) complex 4 were comparable with cisplatin, while complexes 5 and 6 had only a minor effect on activation of these caspases.     

Synthesis of α-fluoro-α,β-unsaturated esters monitored by 1D and 2D benchtop NMR spectroscopy

For optimization and control of pharmaceutically and industrially important reactions, chemical information is required in real time. Instrument size, handling, and operation costs are important criteria to be considered when choosing a suitable analytical method apart from sensitivity and resolution. This present study explores the use of a robust and compact nuclear magnetic resonance (NMR) spectrometer to monitor the stereo-selective formation of α-fluoro-α,β-unsaturated esters from α-fluoro-β-keto esters via deprotonation and deacylation in real time. These compounds are precursors of various pharmaceutically active substances. The real-time study revealed the deprotonation and deacylation steps of the reaction. The reaction was studied at temperatures ranging from 293 to 333 K by interleaved one-dimensional ¹H and ¹⁹F and two-dimensional ¹H–¹H COSY experiments. The kinetic rate constants were evaluated using a pseudo first-order kinetic model. The activation energies for the deprotonation and deacylation steps were determined to 28 ± 2 and 63.5 ± 8 kJ/mol, respectively. This showed that the deprotonation step is fast compared with the deacylation step and that the deacylation step determines the rate of the overall reaction. The reaction was repeated three times at 293 K to monitor the repeatability and stability of the system. The compact NMR spectrometer provided detailed information on the mechanism and kinetics of the reaction, which is essential for optimizing the synthetic routes for stepwise syntheses of pharmaceutically active substances.     

Dynamics of laser-induced cavitation bubble during expansion over sharp-edge geometry submerged in liquid – an inside view by diffuse illumination

Laser ablation in liquids is growing in popularity for various applications including nanoparticle production, breakdown spectroscopy, and surface functionalization. When laser pulse ablates the solid target submerged in liquid, a cavitation bubble develops. In case of “finite” geometries of ablated solids, liquid dynamical phenomena can occur inside the bubble when the bubble overflows the surface edge. To observe this dynamics, we use diffuse illumination of a flashlamp in combination with a high-speed videography by exposure times down to 250 ns. The developed theoretical modelling and its comparison with the experimental observations clearly prove that this approach widens the observable area inside the bubble. We thereby use it to study the dynamics of laser-induced cavitation bubble during its expansion over a sharp-edge (“cliff-like” 90°) geometry submerged in water, ethanol, and polyethylene glycol 300. The samples are 17 mm wide stainless steel plates with thickness in the range of 0.025–2 mm. Bubbles are induced on the samples by 1064-nm laser pulses with pulse durations of 7–60 ns and pulse energies of 10–55 mJ. We observe formation of a fixed-type secondary cavity behind the edge where low-pressure area develops due to bubble-driven flow of the liquid. This occurs when the velocity of liquid overflow exceeds ~20 m s⁻¹. A re-entrant liquid injection with up to ~40 m s⁻¹ velocity may occur inside the bubble when the bubble overflows the edge of the sample. Formation and characteristics of the jet evidently depend on the relation between the breakdown-edge offset and the bubble energy, as well as the properties of the surrounding liquid. Higher viscosity of the liquid prevents the generation of the jet.     

An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.